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ABSTRACT 

Over the last few years, the Sloan Digital Sky Survey (SDSS) has discovered several 
hundred quasars with redshift between 4.0 and 6.4. Including the effects of magnification 
bias, one expects a priori that an appreciable fraction of these objects are gravitationally 
lensed. We have used the Advanced Camera for Surveys on the Hubble Space Telescope 
to carry out a snapshot imaging survey of high-redshift SDSS quasars to search for 
gravitationally split lenses. This paper, the first in a series reporting the results of 
the survey, describes snapshot observations of four quasars at z = 5.74,5.82,5.99 and 
6.30, respectively. We find that none of these objects has a lensed companion within 
5 magnitudes with a separation larger than 0.3 arcseconds; within 2.5 magnitudes, we 
can rule out companions within 0.1 arcseconds. Based on the non-detection of strong 
lensing in these four systems, we constrain the z ~ 6 luminosity function to a slope of 
P > —4.63 {3a), assuming a break in the quasar luminosity function at MJ^^q = —24.0. 
We discuss the implications of this constraint on the ionizing background due to quasars 
in the early universe. Given that these quasars are not highly magnified, estimates of 
the masses of their central engines by the Eddington argument must be taken seriously, 
possibly challenging models of black hole formation. 

Subject headings: gravitational lensing — early universe — quasars: general — quasars: 

individual (SDSSp J104433.04-012502.2, SDSSp J083643.85+005453.3, SDSSp J130608.26+035626.3, 

SDSSp J103027.10+052455.0) — galaxies: luminosity function 
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1. Introduction 

Before 1995, only 24 quasars with redshifts larger than four had been published; most had been 
found in multicolor or grism surveys for high-redshift quasars (e.g., Waxren et al. 1987; Schmidt, 
Schneider, & Gunn 1987). In the past eight years the number of known z > 4 quasars has increased 
to more than 400, largely due to the discoveries from the Sloan Digital Sky Survey (SDSS; York 
et al. 2000) and the Digitized Palomar Sky Survey (Kennefick, Djorgovski, & Meylan 1996). The 
SDSS to date has published the discovery of over 200 quasars with redshifts greater than 3.6 (Fan 
et al. 2001b; Anderson et al. 2001; Schneider et al. 2001; Schneider et al. 2003). In addition. 
Fan et al. (2000), Fan et al. (2001a), and Fan et al. (2003) have discovered the seven highest- 
redshift quasars known, all with z > 5.7 and four of which are the subject of this paper. This large 
sample provides us with new opportunities to study high-rcdshift quasars; it has already been used 
to determine the luminosity function of quasars and its evolution with rcdshift from z = 3.6 to 
z = 6.0 (Fan et al. 2001c, 2003). The luminosity function is well-fit with a power-law, while the 
comoving number density of quasars drops exponentially with redshift over this range. One of the 
goals of this paper is to determine the extent to which the observed luminosity function might be 
biased by gravitational lensing. 

Understanding the (intrinsic) quasar luminosity function and evolution thereof at high redshifts 
is important for several reasons. We observe that supermassive, quiescent black holes are ubiquitous 
at low redshift (Magorrian et al. 1998); comparison of the present-day mass function with the high- 
redshift luminosity function can constrain models for the duty cycle of black holes and their feeding 
mechanisms. There is increasing evidence for a population of optically faint AGN at high redshift, 
which appear in deep X-ray images with the Chandra satellite (Mushotzky et al. 2000; Brandt et al. 
2001; Alexander et al. 2001); understanding their nature requires that we understand the nature 
of the optically luminous objects as well. Finally, the optical (rest-frame UV) luminosity function 
puts important constraints on the ultraviolet ionizing background (Fan et al. 2001c), which is 
especially important now that we appear to be observing the first hints that the universe became 
completely reionized around redshift z = 6 (e.g., Becker et al. 2001; Djorgovski et al. 2001; White 
et al. 2003). 

The continuum luminosities of the z > 5.7 quasars arc all very high, with absolute magnitudes 
at rest frame 1450 A between —26.8 and —28 for the WMAP cosmology that we use throughout 
this paper (a fiat A CDM model with = 0.29 and Hq = 72 km s~^ Mpc~^; Spergel et al. 
2003) . This yields a lower limit to the masses of their black holes of order 4 x 10^ Mq by the 
Eddington argument (Fan et al. 2000). It is a challenge to explain in standard cosmologies how 
such massive black holes could form at a cosmic epoch less than a billion years after the Big Bang 
(e.g., Efstathiou & Rees 1988; Turner 1991; Haiman & Loeb 2001). 

All these results assume that the apparent magnitudes that we measure for these quasars are 
intrinsic to the quasar themselves, and represent isotropic fluxes. However, gravitational lensing 
can systematically brighten the images of quasars by a significant amount. This is a particu- 
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larly important effect for flux- limited surveys like SDSS, which are sensitive to the most luminous 

quasars at any given redshift (e.g., Maoz & Rix 1993). The resulting magnification bias (e.g., Pci 
1995, and references therein) can systematically change the quasar luminosity function, causing an 
overestimation of the black hole masses powering these objects. 

The expected fraction of multiply imaged quasars at a given redshift depends both on the 
cosmological model (which sets the distribution and profiles of halos as a function of redshift, as 
well as the dependence of cosmological pathlength on redshift), and on the luminosity function 
of quasars, which determines the effect of magnification bias. In the WMAP cosmology (using a 
halo distribution taken from the large N-body simulations of Jenkins et al. 2001), the fraction 
of random lines of sight aX z = 4 that produce multiple images at all splitting angles (see the 
discussion of lensing models in § 3 below) is of order 0.2%; this fraction rises to 0.4% at z = 6. 
However, magnification bias greatly increases the lensing rate in a flux-limited sample. At z = 4, 
using the canonical cosmology, and a luminosity function $(-L) oc (Fan et al. 2001b) with 

a break at Mj*45Q = —24.0 (see § 3 below for further discussion), we expect about 1% of quasars 
in our sample to be lensed (specifically, multiply imaged with magnification of at least two) . This 
fraction increases to 8% if we assume a steeper slope of /? = 3.5, and to ~ 2% if we keep the slope 
of /3 = 2.5, but assume a break in $ at a 10 times lower luminosity {M*^^q = —21.5). At z = 6, 
magnification bias can result in lensing probabilities approaching 100% if the luminosity function is 
particularly steep (Comerford, Haiman, &; Schaye 2002; Wyithe Sz Loeb 2002a, 2002b). However, 
the detailed shape of $ is not well-measured, especially at higher redshifts and luminosities. Such 
quasars represent the highest peaks in the density field at high redshift, where the mass function 
of massive dark halos — in which bright quasars likely reside — is expected to be very steep (Fan 
et al. 2001c). 

The median expected splitting of gravitationally lensed quasars is predicted to be somewhat less 
than 1" (Turner, Ostriker, & Gott 1984), with flux ratios as large as 20:1; this is in broad agreement 
with the observed splitting separation found by the Cosmic Lens All Sky Survey (CLASS; Browne 
et al. 2003), but see the discussion in §3 below on CLASS implications for the details of the lens 
profiles. The SDSS images themselves have PSF widths of order 1" — 1''5 (Abazajian et al. 2003), 
making them inadequate to look for all but the very rare wide-separation (A^ > l'.'5) lenses (Pindor 
et al. 2003). We are thus carrying out a snapshot imaging survey of 250 SDSS-discovcrcd quasars 
with z > 4.0 using the Advanced Camera for Surveys (ACS; Ford et al. 2003) on the Hubble Space 
Telescope (HST), which can easily resolve pairs separated by as little as O'.'l with flux ratios as large 
as 10:1. Our observations should therefore be sensitive to essentially all expected lenses (Turner 
et al. 1984; Hinshaw &; Krauss 1987). Fan et al. (2003) presented some preliminary results based 
on some of these (and other, non-HST) data; here we give a full analysis of the HST observations 
for all four of the z > 5.7 quasars in our sample. 

Throughout this paper, we adopt the WMAP cosmology mentioned above, with an rms mass 
fluctuation within a sphere of radius 8 Mpc of dg = 0.9, and power-law index n = 0.99 for 
the power spectrum of density fluctuations (Spergel et al. 2003). We also adopt the cosmological 
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transfer function from Eisenstein & Hu (1999). Conversions between Mb and M1450 assume Mb = 
Mu5o — 0.48 (Schmidt, Schneider, & Gunn 1995) with spectral index ai, = —0.5 {fi, oc i^"""). 

2. The Data 
2.1. Observations 

The full sample of z > 4.0 quasars included in our HST snapshot program will be presented 
in future papers in this series. In this paper, we present results of imaging the four quasars with 

z > 5.7 discovered by Fan ct al. (2000, 2001a); sec Tabic 1. These four objects were selected from 
the SDSS imaging database based on their extremely red i — z color, absence of detection in u, g, 
and r, and relatively blue z — J color (which distinguishes the z > 5.7 quasars from brown dwarfs). 

The SDSS images of these quasars have an image scale of 0'.'396pixel~^ and z-band seeing 
FWHM between l'.'4 and I'.'S. In order to test the ability of the SDSS photometric pipeline to 
identify quasar pairs, we created simulated SDSS images of pairs of point sources (see Pindor et al. 
2003 for details) in these observing conditions and at the appropriate signal-to-noise ratio. We find 
that the SDSS star-galaxy separator would have identified pairs of point sources as resolved objects 
for image separations greater than ~ 1" and for flux ratios less than ~ 5:1. Thus the SDSS data arc 
suitable for exploring only a fraction of the parameter space of lensing (in terms of separations and 
flux ratios) that is of interest; higher resolution HST observations were needed to probe the rest of 
that parameter space. We did not restrict our sample to point sources since the SDSS star-galaxy 
separation algorithm (see Scranton et al. 2002) breaks down at these faint magnitudes. The only 
object of the four that was flagged as an extended source was SDSSp J104433.04-012502.2. 

We used the High Resolution Camera (HRC) on ACS to image the four z > 5.7 quasars 
presented herein. For these four objects, the Lya forest lies in the SDSS z-band; we therefore 
observed them in the SDSS 2;-band, F850LP, despite its lower sensitivity. The Wide-Field Camera 
on ACS has higher sensitivity in i and z than does the HRC, but has has substantially higher 
overhead and is mildly undersampled. The exposure times were 1200 seconds for each object, 600 
seconds in each of two exposures to help in cosmic ray rejection. We did not dither the images, and 
our cosmic ray rejection is far from perfect (§ 2.2); however, our analysis indicates that imperfect 
cosmic ray removal is not critical for our search for gravitational lenses. 
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2.2. Data Processing 

The raw images were calibrated by the CALACS package in IRAF^ as part of on-the-fly- 
reprocessing (OTFR) at the time of download. For the most part, we find these initial reductions 
to be sufficient for our purposes. The images that we present are the "cosmic ray rejected" (CRJ) 
images that arc output by the OTFR algorithms at STScI. These files were re-requested a number of 
weeks after the observations were completed to ensure that we are using the latest "reference" files 
and that the OTFR reprocessing is current. The CRJ files have all been reduced in the standard 
manner, including having been overscan-, bias- and dark-corrected, flat-fielded and photometrically 
calibrated, in addition to having bad pixels masked and cosmic rays removed (see the ACS manual^ 
for more details). The plate scale on the HRC of the ACS is not the same in the x and y directions 
and the sky axes are not perpendicular on the CCD; these CRJ images have not been corrected for 
this distortion. A 5" x 5" cutout of each of the ACS images is given in Figures 1-4. 

With the exception of cosmic ray rejection and hot pixel identification, the images output 
by the OTFR algorithm are acceptable for our purposes. The CALACS rejects pixels that are 
significantly different in the two images. This process successfully rejects the majority of the 
cosmic rays in our images. However, many of the cosmic rays occupy multiple pixels, and we have 
only two exposures, meaning that some cosmic rays remain. In addition, hot pixels can remain in 
the corrected images if they have not been properly identified in the masks that are constructed 
for each period of observations between the annealing procedure that helps to repair hot pixels. 
Fortunately, the morphologies of the residual cosmic rays and hot pixels are considerably different 
from that of the ACS PSF (the PSF is spread out over many tens of pixels), so spurious sources 
due to cosmic rays and most other defects are easy to recognize. 

2.3. Subtracting the Point Source 

The point spread function is narrow enough that any lens with a separation greater than ~ ^'.2 
will be obvious by visual inspection (see § 2.5). To search for faint secondary images at separations 
smaller than this requires that we fit and subtract off a model for the point spread function of each 
image. We have done this using the v6.0 (pre-launch of ACS) version of the Tiny Tim software 
(Krist 1995)^, which produces a model PSF for the instruments on HST given the object's SED and 
position in the focal plane, the filter curve, and knowledge of the optics of the instrument. Such 
model PSFs are necessary because it is not always possible to do PSF subtraction using an empirical 
template (another point source in the field-of-view). That is, given the positional dependence of 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 

*http: / / www.stsci.edu /hst / acs / documents /handbooks/cyclel2/cover.html 

^http: / /www. stsci.edu/software /tinytim/ 
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the PSF across the field-of-view, the Tiny Tim models are often better than using an empirical 
template. For our ACS/HRC images, the model correctly takes into account the different scales 
and non-orthogonality of the x and y axes. 

We fit the PSF model to each CRJ image, allowing the location on the CCD and normalization 
to vary, and using sine interpolation when shifting the model PSF by fractional pixels. Inputs to 
Tiny Tim were the camera used (ACS/HRC), the location on the CCD of each quasar to the nearest 
pixel (as measured by imexam in IRAF) and a Keck spectrum (Becker et al. 2001) of each quasar 
(for the wavelength dependence of the PSF). We minimized the sum of the square of residuals 
(since the formal errors per pixel over the image are close to constant). The results are shown in 
Figures 1-4. Each of the four images show the familiar first Airy ring; on a very hard stretch, 
the second Airy ring is faintly visible. Note that HST PSFs are not symmetric as a result of the 
three primary mirror supports, the signature of which can be scon in the PSF subtracted images. 
Given inevitably somewhat imperfect PSF subtraction and the expected asymmetry of the residual 
images, in none of the four cases is there an obvious secondary image; none of them are obviously 
gravitationally lensed. 

2.4. Notes on Individual Objects 

24.1. SDSSp J104433.04-012502.2; z = 5.74 

The ACS/HRC image of this quasar is shown in the left panel of Figure 1; the right panel 
shows the residuals after PSF subtraction. The SDSS z-band magnitude of this quasar is 19.3; 
the HST 2-band (ST) magnitude is 19.91, see also Table 1. The offset between the SDSS and 
ST z-band magnitude is systematic (sec the other objects below) with the ST magnitudes always 
being fainter. This is probably due to the differences in the filter curves (the SDSS curve having 
less red response due to atmospheric features) and the differences between a Vega-based (ST) and 
an AB-based (SDSS) photometric system given that Vega has an AB z-band magnitude of 0.57 
(Fukugita et al. 1996). We find no offset between the SDSS and ST z-band magnitudes (in the 
HST images of our lower redshift quasars), where the filter curves appear to be much more similar. 

Fan et al. (2000) describe a K' image of this object taken with the Keck 10m telescope in 
0'.'375 seeing, which was consistent with a point source. That is consistent with the HST image; 
there is no apparent splitting in this object. 

Shioya et al. (2002) found a faint galaxy {ruB ^ 25) 1''9 south-west of this quasar. Our 20- 
minute exposure shows this object faintly after rebinning. Shioya et al. (2002) suggest that this 
galaxy could amplify the flux of the quasar by up to a factor of two. The lack of a secondary image 
of the quasar restricts any possible magnification to be no larger than this — in the absence of 
microlensing and external shear (Wyithe & Loeb 2002a). 
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2.4-2. SDSSp J083643.85+005453.3; z = 5.82 

The ACS/HRC image of this quasar and residuals after PSF subtraction are shown in Figure 2. 
The SDSS z-band magnitude of this quasar is 18.79; the HST z-band (ST) magnitude is 19.49, see 
also Table 1. One other point source and two other extended sources are in the ACS/HRC field 
of this quasar, but all are more than 5" away (out of the ficld-of-vicw of Fig. 2), and are therefore 
unlikely to lens (or significantly magnify) the quasar. The black dots in the field-of-view of Figure 2 
are residual cosmic rays that were not properly cleaned by the OTFR reductions. 

2.4.3. SDSSp J130608.26+035626.3; z = 5.99 

The ACS/HRC image of this quasar is shown in Figure 3. The SDSS z-band magnitude of 
this quasar is 19.46; the HST z-band (ST) magnitude is 20.23, see also Table 1. There is some 
evidence for a few very faint extended sources in the ACS /HRC field (off the image in Fig. 3) , but 
as with Figure 2, all of the black dots in Figure 3 are residual cosmic rays. Yamada et al. (2003) 
have taken deep images to i ^ 25.5 in moderate seeing with Subaru in this field, and have found 
no foreground galaxies close enough to the line of sight to cause gravitational lensing. 

2.4.4. SDSSp J103027.10+052455.0; z = 6.30 

The ACS/HRC image of this quasar is shown in Figure 4. The SDSS z-band magnitude of 
this quasar is 20.02; the HST z-hand (ST) magnitude is 20.54, see also Table 1. There is one other 
relatively bright point source in the field, nearly due west of the target quasar, but it is at the edge 
of the field (off the image in Fig. 4) and is not of interest in terms of lensing. As in the previous 
object, a deep image by Yamada et al. (2003) shows no obvious lensing galaxies in the foreground. 
However, Petric et al. (2003) have taken deep VLA images at 20 cm in this field, and found a 
statistically significant overdensity of four radio sources within l' of the quasar, which may indicate 
the possibility of a magnifying cluster along the line of sight. 

2.5. Looking for Multiple Images 

Our ACS images have an image scale of 0'.'025pixel~^, and we can visually rule out any 
multiple imaging at scales larger than about 0'.'2. To make this more quantitative, we have carried 
out simulations of quasar pairs at various separations and flux ratios, by adding a scaled-down and 
shifted version of our quasar images to itself. This naive procedure produces a slightly incorrect 

sky level, but the inaccuracy is negligible for our present purposes. Figure 5 shows simulated 
observations of SDSS J0836+0054 with an additional point source 2.5 and 5.0 mag fainter than 
the quasar and separated from the quasar by 0'.'3. The secondary image is obvious in both cases 
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even without PSF subtraction. Figure 6 shows the result of a (/.'l offset. With a flux ratio of 10:1, 

the secondary object is visible as an enhancement in the first Airy ring of the primary object; it 
becomes clear upon PSF subtraction (Figure 7). With a flux ratio of 100:1, we cannot discern a 
pair with O'.'l separation even after subtraction of the PSF; however, even if we were able to do so, 
it would not significantly improve the constraints that we derive below since most split images will 
have larger image separations. 



3. Constraining the Slope of the Quasar Luminosity Function 

The shape of the high-redshift quasar luminosity function is not tightly constrained because 
of the limited size and luminosity range of the available statistical samples. Schmidt et al. (1995), 
using a set of 90 quasars with 2.7 < z < 4.8, found a power-law luminosity function slope (/3) of 
~ —2. Fan et al. (2001c) measured (3 ~ —2.5 it 0.3 from 39 quasars in the redshift interval 3.6 to 
5.0. These slopes apply to the high-luminosity (M < —26) part of the luminosity function. At these 
high redshifts, direct observations currently yield only poor constraints on the slope of the quasar 
luminosity function (LF) at z ~ 5, and essentially no constraints at z 6, where the luminosities 
of the known quasars all lie within 0.5 magnitudes of the detection threshold of SDSS. From the 
number counts of z ~ 6 quasars. Fan et al. (2001c) place a 3-cr lower limit on the slope of /? = —4.6. 
Gravitational Icnsing tends to flatten the intrinsic luminosity function, producing an apparent slope 
of /3 = — 3 in the limit of an arbitrarily steep intrinsic LF down to a fixed low-luminosity cutoff 
(e.g., Schneider, Ehlers, & Falco 1992). As a result, the relatively shallow apparent slope at z ~ 5 
does not necessarily reflect the slope of the underlying LF, which could be much steeper, if most 
quasars were lensed and magnified. 

One of the primary goals of this work is to constrain the slope of the quasar LF at high redshift, 
using the lack of evidence for strong lensing among the four quasars in the HST images (Comerford 
et al. 2002). In line with previous work, we describe the intrinsic (not necessarily observed) quasar 
luminosity function as a broken power law (e.g., Boyle, Shanks, &; Peterson 1988; Pei 1995): 

- {L/L^Y^ + {L/L^f^- 

The LF is described by four parameters: the normalization the faint-end slope /?/, the bright-end 
slope f3h, and the characteristic luminosity L* at which the LF steepens. The lensing probability 
is very sensitive to the last two parameters, (3^ and L*, while it is only weakly dependent on the 
faint-end slope below the break, and it is strictly independent of the normalization Here we 
apply the lensing model from Comerford et al. (2002), in which lenses are associated with dark 
matter halos, to compute the total lensing probability, including the effect of magnification bias. 
In this model, the abundance of halos as a function of potential well depth is adopted from the 
simulations of Jenkins et al. (2001). While it would be simplest to assume that the halos have a 
common density profile, such as the singular isothermal sphere (SIS), recent work has suggested 
that the typical density profile of massive lenses is shallower than that of the SIS (Porciani & 
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Madau 2000; Kochanek & White 2001; Li & Ostriker 2002; Huterer &; Ma 2003). The observed 
dearth of large-separation lenses cannot be explained if all lenses have the same density profile. 
A prescription that fits the observed steep splitting angle distribution is to assume that all halos 
below M 10^^ Mq have SIS profiles (adopting a standard conversion between circular velocity 
and halo mass), while all halos above this mass follow the dark matter density profile suggested by 
Navarro, Prenk, &; White (1997, hereafter NFW). NFW profiles are much less efficient lenses than 
are SIS profiles. This prescription is physically motivated in Kochanek & White (2001) based on 
the efficiency of gas cooling, and we adopt it in this paper. In our calculations, this is essentially 
equivalent to ignoring all lenses above a halo mass 

of lO^^M^^ as the massive halos do not contribute 

to lensing at small separations. 

As an example, we first fix the characteristic luminosity ( "knee" ) of the quasar LF at a fiducial 
value of M^45Q = —24.0, and derive constraints on the bright-end slope, Ph- The choice for this 
value is somewhat arbitrary: the highest redshift where the quasar LF is measured to faint enough 
magnitudes to define a reliable break is at z 2.15 (Boyle et al. 2000), where M^^^q ^ —25.6 (Mg 
—26.1). In pure luminosity evolution models, the break at z ~ 5 would then be at around M^^^q = 
-24.0 (e.g. Madau, Haardt, & Rccs 1999; Wyithc & Loeb 2002a). Having fixed the value of the 
break, we next vary the slope of the LF, and require that the expected probability for no lensing 
among the four highest redshift quasars is P=0.0026 (= 3a) or greater. Assuming that we can rule 
out a second image with a flux ratio of 100 down to a splitting angle of O'.'S, we find the constraint 
on the slope /J^ > —4.63. We find that resolving high-contrast lensing configurations, with a flux 
ratio > 10, or resolving angular separations smaller than 0'.'3 does not significantly improve this 
constraint, since in any model, the fraction of lenses with separations less than O'.'S is small. For 
example, ruling out flux ratios of 10 down to O'.'l would only improve the constraint to Ph > —4.58. 

It is more informative to illustrate joint constraints on the parameters {M*^^q, (3fi)- In Figure 8, 
we show contours for the likelihood at P = 0.32, 0.05, 0.01, and 0.001 that none of the four z > 5.7 
quasars are lensed (again, to a separation of 0.3 arcseconds, and to flux ratios of 100). The figure 
shows, as in Fan et al. (2003), that the derived limit on the slope Ph depends strongly on the 
assumed break. Overall, the constraint shown in this figure is similar to that derived from the 
sample of seven high-redshift quasars in Fan et al. (2003), although the present constraints are 
somewhat weaker, due to the smaller number of sources. The recent results of Wyithc (2003) are 
consistent with ours under the assumption of no lensing. As shown in Wyithc (2003), a stronger 
constraint is available if two of the known z > 6 quasars are indeed magnified, as suggested by 
Shioya et al. (2002) and White et al. (2003). 



4. Discussion 

The lack of lenses among the four high-redshift quasars allows us to place constraints on the 
quasar LF at high redshift, shown in Figure 8. It is useful to keep in mind that only models with 
extreme magnification bias, predicting significant lensing probabilities, can be constrained by the 
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kcfc of lensing events. For example, for the range of LF slopes considered by Wyithe & Loeb (2002a), 
f3h = 2.58 — 3.43, the lensing fraction for the z ~ 6 SDSS sample would be 7-30%. The chance 
of finding no lenses among the four quasars for their steepest slope of 3.43 is (1 — 0.3)^ = 0.24, 
which translates to a significance of only ~ la for ruling out this slope (consistent with our results 
shown in Figure 8 at the break of M1450 = —24.0). 

Accordingly, as Figure 8 shows, the constraints obtained here still allow very steep LF slopes, 
especially if we assume that the break in the LF lies at relatively bright magnitudes. Nevertheless, 
these constraints are fairly direct, and are among the most stringent limits we currently have on 
the shape of the high-redshift quasar LF. We expect that the full z > A sample will further tighten 
the constraints (or yield a measurement should a lens be found). Table 2 shows that the lensing 
probability (for a slope of f3 = —3.5 and break M1450 = —24.0) decreases from 8.2% at 2; = 6 to 
1.2% at z = 4. This implies that seeing no lensing with a single z = 6 quasar provides a constraint 
that is equivalent to about seven z = A quasars, from the equation (1 — 0.012)^ = (1 — 0.082). 
The lensing probability is a steep function of l3h and M*, but since our full sample will include 
~ 250/4 « 60 times more sources than utilized here, we expect that we will be able to significantly 
tighten the constraints in the {Ph,M*) plane with this sample. 

Given our assumed shape of the quasar LF, the constraints obtained in the (/3/i,M*) plane 
translate directly to the global ionizing emissivity of the z ~ 6 quasar population, uei, = J L^{L)dL. 
The emissivity is of great interest, since a related quantity — the mean background flux — can 
be constrained independently from the level of resonant neutral hydrogen absorption detected in 
high-resolution spectra of the SDSS quasars. Specifically, recent upper limits on this background 
were obtained from the Gunn-Peterson (1965) trough of the z = 6.30 quasar (e.g., McDonald & 
Miralda-Escude 2001; Cen k McDonald 2002; Lidz et al. 2002; Fan et al. 2002), requiring that 
r < 10~^^ ionizations per hydrogen atom per second at z = 6. The conversion to the ionizing 
emissivity i>€i, is 

Tvan , (1 + zY ueu 1 , . 

i = -7— Amfp — s , (2) 

/ip Ait i/Yi 

where we have assumed ej, oc v~^, an = 6.3 x 10~^^ cm""^ is the hydrogen ionization cross section 
at the Lyman limit (and we assumed <th{i^) oc above the threshold), hp is the Planck constant, 
= 3.29 X 10^'''Hz is the frequency of a Lyman limit photon, and A^fp is the mean free path of 
ionizing photons (see the Appendix in Comerford et al. 2002). The last quantity unfortunately 
has a very large uncertainty. Nevertheless, for a crude estimate, we here assume that the effective 
mean free path corresponds to a redshift interval of Az = 0.17 at 2; = 3 (see Haardt Sz Madau 1996; 
Steidel, Pettini, Sz Adelberger 2001), and scales with redshift as Amfp oc (1 + z)~^ towards higher 
redshift (Cen &; McDonald 2002). In this case, the model with M^^^g = -24.0 and I3h = -4.6 
corresponds to F 10~^^sec~^. In other words, the constraint we derive here from the absence 
of lensing is approximately at the same level as that implied by the upper limit on the ionizing 
background. Note that the constraints we obtain apply to the total emissivity of quasars alone, but 
the upper limits on the background from the Gunn-Peterson trough constrains the total emissivity 
(including galaxies). Future improvements of the lensing constraint can yield tighter limits on the 
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ionizing emissivity of the quasar population than hitherto possible — subject to the uncertainties 
in (1) the mean evolution of the mean free path with redshift, (2) the faint -end slope which 
does affect the total emissivity ve,^, and (3) the escape fraction of ionizing photons at 2; ~ 6. 

The fact that none of the four high redshift quasars appears to be strongly magnified also 
implies that the masses inferred for the central black holes from the Eddington limit argument are 
real. If a quasar is strongly beamed towards us, the bolometric luminosity would be overestimated, 
but Willott, McLure, Sz Jarvis (2003), following Haiman & Cen (2002) argued that emission line 
equivalent widths being normal mitigates the beaming hypothesis. The inferred black hole masses 
in these four quasars are very high, 2 — 4 x 10^ M0; how black holes can grow so large so early in 
the Universe requires explanation (Efstathiou & Rees 1988; Turner 1991). Haiman &; Loeb (2001) 
considered the presence of the ~ 4 x 10^ Mq black hole for the z = 5.75 quasar in simple models for 
black hole growth. They found that it could be marginally accommodated in simple models of black 
hole growth by mergers and accretion in an LCDM universe. These considerations depend strongly 
on the redshift, and on the typical radiative efficiency of black holes. The highest redshift quasars 
therefore represent a bigger challenge to models, especially if their radiative efficiency during their 
growth was as high as 20%, the typical value inferred recently for the brightest quasars at lower 
redshifts (Yu &: Tremaine 2002). This is twice the fiducial value of ~ 10% adopted in Haiman &: 
Loeb (2001) and most other previous work on black hole growth, doubling the expected growth 
timescale. 

Finally, the lack of strong lensing appears consistent with the large apparent sizes of ionized 
regions around the z 6 quasars. Haiman & Cen (2002) argue that the z = 6.30 source cannot 
be magnified by more than a factor of a few, since an intrinsically faint source would not be able 
to ionize a neutral IGM out to the observed distance away from the source of ^ 30 (comoving) 
Mpc — the spectrum shows transmission at wavelengths corresponding to this distance. Similar 
conclusions were drawn for the z = 6.42 quasar (Fan et al. 2003) recently by White et al. (2003). 

5. Conclusions and Future Work 

We have obtained high resolution HST images of the four highest redshift quasars known prior 
to 2003 to look for the signature of gravitational lensing. We have found no evidence of multiple 
images, significantly limiting the amount by which these quasars can be magnified by foreground 
mass concentrations (in the absence of microlensing). Thus the masses of these quasars as derived 
from the Eddington argument cannot be significantly inflated by magnification. 

The lack of any strong lenses puts a 3cr constraint on the bright end slope of the z > 5 luminosity 
function of Ph > —4.63. This limit on the slope of the luminosity function is also consistent with 
previous limits on the ionizing emissivity of 2; ~ 6 quasars of F ~ 10~^^ sec~^ ionizations per 
hydrogen atom per second based on the Gunn-Peterson trough in the highest redshift quasar 
studied herein. 
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The next paper in the series will present an analysis of a much larger sample of z > 4 quasar 
observations with HST and will discuss the constraints that we can place upon the quasar luminosity 
function from that larger (but lower redshift) sample. Furthermore, there are many more z > 4 (and 
even a few z > 5.7) quasars known now than were known at the time that these i/ST observations 
were approved. Higher resolution images should be taken of as many of these as possible; much 
of the parameter space can even be covered by ground-based observations since HST resolution 
splittings are only needed to discern the smallest expected splittings. 
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Fig. 1.— HST/ ACS /RRC F850LP images of SDSS J1044-0125 {z = 5.74). The scale is 5" x 5" in 

each of the panels. The left hand panel is the "cosmic ray rejected" (CRJ) output of CALACS. The 
right hand panel is the same CRJ image after subtraction of the Tiny Tim PSF (v6.0). Note that 
the 25 faint galaxy that Shioya et al. (2002) detected 1'.'9 to the south-west of the quasar 

position is not visible at the depth of this image. 
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Fig. 2.— As in Figure 1, for SDSS J0836+0054 {z = 5.82) 
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Fig. 4.— As in Figure 1, for SDSS J1030+0524 {z = 6.30). 



Fig. 5. — Simulated lenses using the SDSS J0836+0054 images. {Left) Lens simulated by scaling 
down the PSF of SDSS J0836+0054 by a factor of 10 and shifting it down by (X'3. (Right) Lens 
simulated by scaling down the PSF of SDSS J0836+0054 by a factor of 100 and shifting it down 
by 0'.'3. 



Table 1. 



SDSSp J 


Redshift 


SDSS z 


HST z 


Reference 


104433.04-012502.2 


5.74 


19.3 


19.91 


1 


083643.85+005453.3 


5.82*^ 


18.79 


19.49 


2 


130608.26+035626.3 


5.99 


19.46 


20.23 


2 


103027.10+052455.0 


6.30'' 


20.02 


20.54 


2 



Note. — Column (1) lists the J2000 coordinates, column (2) the 
redshift, columns (3) and (4) give SDSS and HST z-h&nd magnitudes 
and column (5) gives the discovery references, which are are (1) Fan 
et al. (2000) and (2) Fan et al. (2001a). 

•^Stern et al. (2003) find z = 5.77 from CIII] emission. 

''Corrected from z = 6.28 using a fit to templates that allow for 
emission line blueshifts (Richards et al. 2002). 
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Fig. 6. — Simulated lenses using the SDSS J0836+0054 images. (Left) Lens simulated by scaling 
down the PSF of SDSS J0836+0054 by a factor of 10 and shifting it down by O'/l. (Right) Lens 
simulated by scaling down the PSF of SDSS J0836+0054 by a factor of 100 and shifting it down 
by O'/l. 
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Fig. 7. — The same as Figure 6, but with PSF subtraction of the brightest component. {Left) Lens 
simulated by scahng down the PSF of SDSS J0836+0054 by a factor of 10 and shifting it down 
by (y.'l. (Right) Lens simulated by scaling down the PSF of SDSS J0836+0054 by a factor of 100 
and shifting it down by O^'l. For the case on the left, we can clearly see the second source, despite 
the PSF residuals (note the triangular pattern due to the three primary mirror supports). The 
second image on the right is not possible to resolve from the primary image without better PSF 
subtraction (if at all). 
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Fig. 8. — Contours of fixed likelihood for no lensing among the four z > 5.7 quasars, shown in the 

two-dimensional parameter space of the slope and break of the z ^ 6 quasar luminosity function. 
The vertical dotted line shows the 3a limit on /3 given an assumed break in the luminosity function 
of Mj*45o = —24.0 {horizontal dotted line). 
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Table 2. 



z 


4 


5 


6 


Mi450 


-25.48 


-26.05 


-26.51 


^1*450 


-24.73 


-24.42 


-24.12 


diff 


0.75 


1.63 


2.39 




0.21% 


0.29% 


0.37% 


Prob. 


1.2% 


3.5% 


8.2% 


Note. 


— The 


rows are 


(1) the 



absolute magnitude at z = 20.2 
mag, (2) L* from the Wyithe & 
Loeb (2002a) quasar LF, (3) the 
difference of the above two, (4) the 
multiple lousing optical depth (i.e., 
fraction of lines of sight subject to 
lensing), and (5) the lensing opti- 
cal depth with bias, assuming slope 
= 3.5. 



